The Astrophysic al Journal, 713:1282-1291,2010 April 20 

Preprint typeset using MPgX style emulateapj v. Q?>I{)1IQ1 



O 
(N 

Oh- 
< 

in 

(N 



o 

is: 



> 

00 
CO 

o 
o 



DOWNFLOWS IN SUNSPOT UMBRAL DOTS 
A. Ortiz\ L.R. Bellot Rubio^, and L. Rouppe van der Voort^ 

^ Institute of Theoretical Astrophysics, University of Oslo, P.O. Box 1029 Blindern, N-0315 Oslo, Norway; ada@astro.uio.no, and 
2 Instituto de Astrofisica de Andalucia (CSIC), Apdo. 3004, 18080 Granada, Spain 
Received 2010 February 13; accepted 2010 March 9; published 2010 April 1 

ABSTRACT 

We study the velocity field of umbral dots (UDs) at a resolution of Q'.'IA. Our analysis is based on full 
Stokes measurements of a pore taken with the Crisp Imaging Spectro-Polarimeter at the Swedish 1 m Solar 
Telescope. We determine the flow velocity at different heights in the photosphere from a bisector analysis 
of the Fe I 630 nm lines. In addition, we use the observed Stokes 2, U, and V profiles to characterize the 
magnetic properties of these structures. We find that most UDs are associated with strong upflows in deep 
photospheric layers. Some of them also show concentrated patches of downflows at their edges, with sizes of 
about velocities of up to 1000 m s~\ and enhanced net circular polarization signals. The downflows 

evolve rapidly and have lifetimes of only a few minutes. These results appear to validate numerical models of 
magnetoconvection in the presence of strong magnetic fields. 
Subject headings: sunspots — Sun: surface magnetism 



1. INTRODUCTION 

At high angular resolution, sunspot umbrae exhibit small 
bright features embedded in a darker, smoothly varying back- 
ground. These features are called umbral dots (UDs) and oc- 
cur in^ssentially all sunspots but also in pores (for a review, 
see lSobotkalll997|) . Frequently a distinction is made between 
central and peripheral UDs based on their location within the 
umbra (e.g. Grossmann-Doerth 1986; Riethmiiller et al. 2008; 
[Sobotka & Jurcak 2009). Central UDs do not show measur- 
able vertical flows and have similar magnetic field inclinations 
as the surrounding umbra. In contrast, upflows and more hor- 
izontal magnetic fields are detected in peripheral UDs. 

The importance of UDs lies in the fact that they may be 
the signature of convection in sunspots. We know that con- 
vective motions must exist becaus e radiation alon e cannot ex- 
plain the brightness of the umbra (iDeinze The ques- 
tion is whether these motions are field free or magnetized. 
Recent simul ations of magnetoconvec tion in a strong back- 
ground field ( Schiissler & Vogler'2006) suggest that the con- 
vective energy transport inside the umbra is achieved in the 
form of narrow plumes of rising hot plasma and strongly re- 
duced magnetic fields. The plumes are associated with bright 
structures that share many similarities with real UDs, includ- 
ing their sizes (about 300 km) and lifetimes (30 minutes). The 
simulated UDs tend to be elliptical in shape, and many have a 
central dark lane^^ 

According to lSchussler & VogleJ (l2006h . the pile-up of ma- 
terial at the top of the plumes increases the gas density and 
moves the r = 1 level toward higher layers, where the temper- 
ature is lower. This produces the dark lanes of the simulated 
UDs. The upflows attain maximum velocities of 3-4 km s~^ 
just below the solar surface. Near r = 1 they turn horizontal 
and move in the direction of the dark lane, until the gas re- 
turns to deeper layers along narrow channels. The downflows 
occur mainly at the endpoints of the dark lanes, but also in 
the immediate surroundings of the UDs. They usually appear 
as tiny patches of concentrated flows with velocities of up to 
1200 ms~^ at z = km. The velocities at r = 1 are smaller 
because the large opacity of the gas hides the layers where the 
flows are stronger (B harti et al. 2010). 

The simulations of ISchiissler & Vogled (l2006|) suggest that 



UDs are a natural consequence of convection in a strong 
magnetic field, which seems to favor the monolithic sunspot 
model over the cluster model. However, the basic physical 
process at work in the simulations, namely overturning con- 
vection, has not yet been confirmed observationally. Upflows 
are known to exist in UDs; what is missing is a clear, unam- 
biguous detection of return downflows. 

The velocity field of UDs has been studied by a number 
of authors. Most of them report negligible velocities in cen- 
tral UDs (e.g., LS chmidt & Balthasar 1994; Rimmele 20Q8|; 
iRiethmuller etaDl2008l: ISobotka & Jurcak 2009) and upflows 
ranging fr om 100 m s~^ to 1000 m s~^ in peripheral U Ds (e.g., 
Socas-Nav arro et al.ir2 004: Rimmele 2004; Riethmiil ler et al.1 
i2008i: ,Sobotka & Jurcakli2009.) . .Bharti et al.. C2007a) claim 
to have measured upflows of 400 ms~^ and downflows of 
300 ms~^ in UDs, using the Universal Birrefringent Filter 
at the Dunn Solar Telescope. This is a remarkable result 
which might have been compromised by the fact that no 
adaptive optics, phase diversity, or any postprocessing tech- 
nique such as speckle reconstruction were available at the 
time of the observation s. In any case, the downflows de- 
scribed by Bharti etaP (12007 a) seem t o be d ifferent from 
those reported by ISchtissler & Vogled (|2006|) in that they 
extend over much larger are as. The sam e happens to the 
downflow patches detected by lBharti et al.l (2009) using Hin- 
ode observations. Interestingly, the spectropolarimetric mea- 
surements of Riethmiiller et al. ( 2008) and Sobotka & Jurcald 
(|2009|) . also from Hinode, do not show downward motions 
despite their better angular resolution (0'/3 versus 0'/6). 

To detect the downflows predicted by the simulations, one 
would like to have extremely high spatial resolution -higher 
than typically achieved by current instruments- and line spec- 
tra to probe the velocity field at different heights in the at- 
mosphere. Until now it has been difficult to fulfill the two 
requirements simultaneously. Here, we present the first spec- 
tropolarimetric observations of UDs approaching a resolution 
of Of'l. This unique data set, acquired at the Swedish 1 m 
Solar Telescope (SST), allows us to investigate the flow field 
and magnetic properties of UDs with unprecedented detail. 
The measurements show localized patches of downflows at 
the edges of UDs, providing firm evidence for magnetocon- 
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vection in sunspot umbrae. 

The paper is organized as follows. The observations and the 
data reduction are described in Section[2l Section[3]deals with 
the morphology, flow field, temporal evolution, and magnetic 
properties of UDs as derived from the observed Stokes spec- 
tra. Finally, in Section |4] we discuss our results and compare 
them with previous works. 

2. OBSERVATIONS AND DATA PROCESSING 

The observations were obtained with th e CRisp Imag- 
ing Sp ectro-Polarimeter (CRISP) at the SST (IScharmer et al.l 
I2003ah on La Palma (Spain). CRISP is based on a dual Fabry 
Perot interferoni eter (FPI) system similar to that described by 
IScharmed (l2006l) . The instrument has been designed to allow 
diffraction-limited observations in the visible. 

CRISP is equipped with three high-speed Sarnoff CCD 
cameras operating at 35 frames per second. Two of them are 
positioned behind the FPI, and the third is located before the 
FPI but after the CRISP pre-filter. This camera is used as an 
anchor channel for image processing and is referred to as the 
wide-band channel. All the cameras take exposures of 17 ms 
and are synchronized by means of an optical chopper. The 
image scale is 0.''071 pixer\ and the field of view (FOV) is 
-71''x71''. The CRISP transmission profile has an FWHM 
of 6.4 pm at 630 nm. 

The polarization modulation is performed by two nematic 
liquid crystal (LC) variable retarders cycling through four dif- 
ferent states. The time needed to change between LC states 
is smaller than the CCD read-out time (10 ms) by use of the 
transient nematic effect. A polarizing beam-splitter in front of 
the narrow-band cameras divides the light into two orthogo- 
nally polarized beams, which results in a significant reduction 
of seeing-induced cross talk (lLiteslll98"7h . 

The data analyzed here were recorded on 2008 June 12 and 
correspond to a pore located at S09 E30 (0 = 31°, /i = 0.85). 
This pore is the naked umbra of a decaying sunspot (AR 
10998) that showed a rudimentary penumbra in previous days. 
CRISP was used to measure the four Stokes profiles of the Fe I 
630.15 and 630.25 nm lines, each sampled at 15 wavelength 
positions in steps of 4.8 pm, from -33.6 to +33.6 pm. In ad- 
dition, a continuum wavelength point was observed. CRISP 
is capable of very fast wavelength tuning (<50 ms) within the 
same spectral region. We recorded 7 exposures per LC state, 
resulting in 28 exposures per wavelength position. The total 
time for a complete wavelength scan of the two lines was 26 s. 

We followed the pore for 23 minutes during which the at- 
mospheric conditions were variable, with moments of excel- 
lent seeing. To improve the spat ial resolution, we used the 
adaptive optics system of the SST (IScharmer et al.l20Q3bh and 
the Multi- Object Multi-Frame B hnd De convolution image 
restoration technique (MOMFBD: lvan No ort et al. 2005). All 
frames in a complete wavelength scan, amounting to 868 im- 
ages per camera, or 2604 images in total (apart from occasion- 
ally dropped frames due to system hick-ups), were included 
in a MOMFBD restoration. The images were divided into 
overlapping 64x64 pixel subfields (4.''5x4.''5) and aU im- 
ages from each subfield were processed as a single MOMFBD 
set. W e refer the reader to van Noort & Roup pe van der Voort 
(l2QQ8h for more details on the MOMFBD processing strate- 
gies on similar multi- wavelength polarimetric scans. 

The restored images were demodulated (pixel by pixel) and 
corrected for telescope polarization using the telescope model 
developed by Selbing (200 5). Details on the CRISP polari- 
metric calibration and the telescope model can be found in 
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Fig. 1 . — Fe 1 630.15 nm intensity profiles and line bisectors of UDs show- 
ing upflows (solid) and downflows (dotted). The corresponding spatial points 
are marked with crosses in Figure [2] The profiles have been normalized to 
their local continua. The bisectors are indicated by asterisks and span the 
range of intensity levels from 0% to 90% . The 80% bisectors give a velocity 
of -1600 m s~^ for the upfiow and +810 m s~^ for the downflow. 

Ivan Noort & Rouppe van derVoortl(l2008h . After demodula- 
tion, the two sets of orthogonal Stokes images recorded by 
the narrow-band cameras were merged. Precise alignment be- 
tween these images is ensured by the pinhole calibration per- 
formed as part of the MOMFBD processing. The final Stokes 
profiles were corrected for residual I to Q,U, and V crosstalk 
by forcing the polarization to be zero in the continuum. An- 
other correction was applied to remove the intensity gradient 
introduced by the CRISP pre filter (FWHM 0.44 nm). 

The noise level in Stokes Q, U, and V amounts to 2.3 x 
10~^, 1.6 X 10~^, and 1.7 x 10~^ of the continuum intensity, re- 
spectiv ely. This is roug hly a factor of 2 higher than the noise 
level of IScharmer et al.l (|2008), who used 4 times more expo- 
sures. We confirmed that the spatial resolution achieved in our 
Stokes maps is close to the diffraction limit of the SST by ana- 
lyzing intensity cuts of the smallest features in the images. We 
conclude that the spatial resolution is basically limited by the 
pixel size, 0.''071, which slightly undersamples the theoretical 
diffraction limit of the telescope (X/D = 0'/l3 at 630 nm). 

2.1. Line parameters 

In this work, we use line bisectors to obtain the velocity 
field of UDs. The bisector is the locus connecting the mid- 
points of horizontal segments between the two wings of the 
intensity profile (see Figure [T] for an example). We deter- 
mined the bisector positions at different intensity levels, from 
0% (line core) to 80% (line wing near the continuum), with 
the help of the bisec_abs_lin .pro routine included in 
the Kiepenheuer-Institut fiir Sonnenphysik IDL library. The 
line core position was computed by means of a parabolic fit 
around the intensity minimum. The calculations were per- 
formed only for Fe I 630.15 nm, since the red wing of Fe I 
630.25 nm is blended with a telluric O2 line. The bisector po- 
sitions were converted into Doppler velocities assuming that 
the pore is at rest. More specifically, we averaged the bisec- 
tors between the 40% and 70% intensity levels within the pore 
and took this as the zero point of the velocity scale. The rms 
fluctuation of the velocity in the pore is 320 m s~\ as com- 
puted from the 80% bisec tor. This is larg er than the values re- 
ported by Franz & Schlichenmaied (12009) . which we attribute 
to the higher spatial resolution of our observations and the fact 
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Fig. 2. — AR 10998 as observed by CRISP on 2008 June 12 at 10:33 UT. (a) Continuum intensity image from the narrow band channel of CRISP. The contours 
outline the edge of the pore and the arrow marks the direction to disk center. The rectangle encloses the region studied in Figure [3] Note the complex network of 
magnetoconvective structures inside the pore. There are UDs with and without dark lanes, (b) LOS velocity map computed from the bisector position at the 80% 
intensity level. The gray scale ranges from -1.5 kms~^ (black) to 1.5 kms~^ (white). Negative velocities indicate blueshifts. 



that we see clear flow structures embedded in the dark umbral 
background. The systematic error of such a velocity calibra- 
tion is not larger than about 150 m s~^ (Rezaei et al. 2006). 

The availability of spectrally resolved line profiles makes 
it possible to probe the velocity field at different heights in 
the photosphere. Roughly speaking, higher intensity levels 
correspond to deeper layers. The continuum is formed at 
around logr = 0, whereas the Fe I 630.15 nm line core is 
formed at logr = -2.9 in the quiet Sun, according to tBalthasad 
(^1985). These values are consiste nt with the d ifferen ce of ap- 
proximately 270 km reported by iGrec et all (l20ld ) between 
the formation heights of the continuum and the line core of 
Fe I 630.15 nm. It is important to realize, however, that 
the intensity observed at a given wavelength is affected by 
the conditions of a large fraction of the atmosphere (e.g., 
ICabrera Solana et aDl2005l and references therein). Thus, any 
formation height actually represents a broad range of optical 
depths. 

To investigate the magnetic properties of UDs we use sim- 
ple line parameters derived from the observed Stokes Q, U, 
and V profiles. These are the mean linear polarization (LP) 
degree. 



LP = 



/^^;[6^(A) + ^^(A)]^/7/(A)dA 



the mean circular polarization (CP) degree, 
/;''[|y(A)|//(A)]dA 



CP: 



and the net circular polarization (NCP), 



NCP= / y(A)dA//c. 



(1) 



(2) 



(3) 



Here, Ar and Ab represent the limits of integration over the 
line, and Ic is the continuum intensity. The NCP accounts 
for the asymmetries in the Stokes V profile, which are due 
to gradie nts in the niagnet i c field and the line-of - sight (LOS) 
velocity (llUing et al.lll97llAuer & Heaslevlll978i) . 

In addition, we have determined the longitudinal 
component of the magnetic field using the center- 
of-gravity ( COG) method (R ees & Semell 119791: 
iLandi degl'Innocenti & Landolfil 120041) . The LOS com- 
ponent of the field is obtained through the relation 



^LOS '■ 



(A+-A_)/2.0 



(4) 



4.67 X 10-13 ^lA^' 

where gi^ is the Lande factor, Aq is the central wavelength of 
the line in A, and X± are the centroids of the right and left 
circularly polarized line components (I±V), calculated as 

/A[/c-(/±y)]dA 



A4 



(5) 



/[/e-(/±y)]dA 

(IUitenbroekll2003h . The COG method is not affected by satu- 
ration effects in the strong field regime. 

The inclination 7 of the magnetic field vector with respect 
to the LOS is determined as 



7 = arccos - 



vT+^-1 



where 



x = 



(6) 



(7) 



This expression is valid for fully split lines 
(Landi degl'Innocenti & Lan dolfil 12004 ). We apply Equa- 
tion (|6]) to six wavelength samples around the peaks of the 
Stokes V profile (three in the red lobe and another three in the 
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Fig. 3. — LOS velocities in an area of the pore with high density of UDs (see the rectangle in Figure[3. Distances are measured in arcsec. The contour lines 
indicate the edge of the pore, (a) Continuum intensity image, (b)-(f) Velocity maps at intensity levels from 0% to 80%. The circles enclose selected UDs and the 
numbers refer to the downflow patches considered in Table [T] 



blue lobe) and take the mean. Once ^los and 7 are known, 
the magnetic field strength B is obtained as 

5 = 5los/|cos7|. (8) 

To reduce the noise, the line parameters LP, CP, and NCP, 
as well as the magnetic field parameters ^los. 7. and have 
been computed separately for the two lines and then averaged. 

3. RESULTS 

3.1. Morphological properties 

Figure [2a) displays a continuum image of the pore as ob- 
served on 2008 June 12 at 10:33 UT. The direction to disk 



center is indicated by the arrow. The pore covers an area 
of around 6" x 9" and shows a complex network of magne- 
toconvective structures (UDs) of different sizes, shapes, and 
brightnesses. In general bigger structures are found in the 
lower half of the pore. Some of the observed UDs present a 
dark lane along their axis. Two of them even resemble coffee 
beans because of their conspicuous dark lanes; one is located 
at (x,y) = (4'', and the other at (x,y) = (5'',5'/5). These two 
specific UDs have a size of O'/l and the smallest dark lane 
is 0'/l4 in widt h (p robably unresolved), in agreement with 
iRimmeld (l2008h and lSobotka & Puschmann (2009). Smaller 
UDs without a dark lane have a typical diameter of 0'/3 and 
minimum sizes of less than Of'2. 
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TABLE 1 
Downflows Observed in Selected UDs 



Patch 


Size 


Distance 


vlos Distance 
(ms-i) 


Vlos 
(m s-i) 


UD size 


/max/^QS 


1 


0'/21 X 0'/2l 


0'/21 


580 


Of' 14 


-1150 


0f'21 x0f'14 


0.71 


2 


0'/21 x0'/21 


Of'28 


680 


0'' 


-620 


Of'70 X Of'35 


0.94 


3 


0f'70 X 0f'21 


Of' 14 


440 


0'' 


-620 


0f'70 X Of'35 


0.94 


4 


Of'28 X Of'28 


0f'21 


990 


0'' 


-380 


Of'28 X Of'35 


1.07 


5 


0f'28x0f'14 


Of'28 


650 


Of'07 


-1420 


Of'28 X Of'28 


1.19 


6 






520 






0f'21 x0f'21 


0.78 


Note. 


— Column 1: 


number of downflow patch in Fi^ 


^urelHf). Column 2: 


size of downflow patch. Column 



3: distance of downflow patch to UD intensity maximum. Column 4: downflow velocity at the 80% intensity level. 
Column 5: distance of upflow patch to UD intensity maximum. Column 6: upflow velocity at the 80% intensity 
level. Column 7: size of UD. Column 8: UD intensity maximum, normalized to the average quiet Sun continuum 
intensity. 



3.2. Velocity field 

The main goal of this work is to obtain the velocity field 
of the UDs present in the pore. Figure [2tb) shows the LOS 
velocities derived from the bisector at the 80% intensity level. 
The gray scale is clipped at ±1.5 kms~^ to emphasize the 
weaker flows. Negative velocities indicate blueshifts. Given 
the heliocentric angle of the observations and the fact that the 
magnetic field is nearly vertical, blueshifts most probably rep- 
resent upflows and redshifts downflows. 

The granulation pattern is clearly visible, while the pore 
shows a complex network of small-scale flow elements. One 
of the striking features of this map is the strong downflows 
seen at the very edge of the pore, near (x^y) = (1" ^9"). They 
will be explored in a forthcoming publication. Most of the 
umbra has small velocities close to zero, but one can detect 
the presence of many upflows and some downflow patches 
inside the pore. These structures are located at the position of 
the UDs observed in Figure Eta), with small deviations that 
we will discuss later. The upflows reach up to -1400 ms~\ 
while the donwflows show a more modest value of around 
400-1000 m s"^ at the 80% intensity level. 

Figure [3] displays the velocity field within the boxed area of 
Figure [21 For convenience, panel (a) shows the corresponding 
continuum image. The circles mark selected UDs, and are 
centered at the position of maximum brightness. The arrows 
indicate the two UDs with dark lanes. Panels (b)-(f) show the 
bisector velocities at five intensity levels, from 0% to 80%. In 
this way, we probe different layers of the photosphere. 

The bisectors close to the line core reveal small Doppler 
shifts. The upflows and donwflows associated with the UDs 
become gradually stronger toward the continuum, such that 
the larger velocities are seen at the 80% intensity level [Fig- 
ure Of)]. Thus, the flows occur preferentially in deep photo- 
spheric layers, and a gradient of velocity with height appears 
to exist. Such a gradient may partly explain the difficulty of 
detecting flows in UDs; it is common to use the line core for 
velocity determinations, but these measurements probe too 
high layers already devoid of flows. 

The majority of UDs show upflows, with maximum veloc- 
ities of around 1400 m s~^ at the 80% intensity level. The 
average size of the upflow patches is around 0.''25, slightly 
smaller than the UDs themselves. It is interesting to note 
that the location of the upflow do not always exactly overlap 
with the position of maximum brightness of the correspond- 
ing UD. Displacements of 0. ''07-0 .''14 (1-2 pixels) are com- 
mon at all intensity levels. This is best seen in Figure Of), 



where the photometric centers of the UDs and the associated 
blue patches tend to be shifted, although the flows always re- 
main inside the circles representing the UDs. 

Downflows are observed in some UDs, but almost exclu- 
sively at the 80% intensity level^ In those deep layers, they 
tend to appear at the edges of UDs exhibiting upflows. How- 
ever, we note that not all UDs with upflows have an associated 
downflow. Part of this may be due to their transient nature 
(Section [33]) . In the case of the UDs showing a dark lane, the 
upflows are cospatial with the dark lane and the downflows 
occur at the endpoints. 

The properties of some selected downflows are presented in 
Table [T] The numbers in the first column identify the corre- 
sponding patches in Figure [3f). All of them are related to an 
adjacent upflow, except for patch 6 which does not seem to 
be associated with any measurable upflow (this patch is em- 
bedded in a large area of redshifts, so it may not be typical). 
The properties of the upflows are given in Table [T] along with 
those of the parent UDs. The average size of the downflow 
patches is about 0."25, but elongated areas of up to 0."7 also 
exist. Since the downflows occur at the periphery of UDs, 
they can be shifted by as much as 0."28 (4 pixels) from the 
position of maximum brightness. There seems to be a pre- 
ferred direction of displacement which coincides with the line 
of symmetry (i.e., the direction to disk center), even though 
not in all cases, since patch 6 does not follow that direction. 
The displacements might simply reflect an enhanced visibility 
of the downflows near the line of symmetry due to projection 
effects. The LOS velocities range from 400 m s~^ to almost 
1000 m s~^ at the 80% intensity level. We do not find any 
relationship between the area of the downflow patch and its 
velocity (patch 4, with the highest measured velocity, is not 
particularly big or small), or between the area and the distance 
to the UD in continuum intensity. 

3.3. Temporal evolution 

We have also studied the temporal evolution of the UDs in- 
side the pore. Their morphological properties change rapidly 
with time. This is demonstrated in Figure [H where six snap- 
shots covering the period from 10:28 UT to 10:49 UT are 
presented. The upper, middle, and lower panels show con- 
tinuum images, velocity maps derived from the 80% bisector, 

^ The extended areas of weak redshifts seen in the line core map coincide 
with the darkest parts of the umbra, so they are not related to UDs. The 
redshifts could be produced either by oscillations or by downward motions of 
the pore as a whole. 
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Fig. 4. — Temporal evolution of UDs in AR 10998. Panels (a)-(f) show six snapshots taken between 10:28 and 10:49 UT on 2008 June 12. Marks have been 
made every arcsecond. The upper panels show continuum images zooming in the same area as Figure [3] The middle and lower panels display velocity maps 
derived from the 80% bisector and NCP maps. The evolution of the UDs indicated by the circles is discussed in the text. 



and NCP maps, respectively. The circles mark prototypical 
UDs whose evolution we have followed in more detail. 

The green circles trace the evolution of a small UD in the 
lower half of the pore. The apparent size of this UD changes 
from a diameter of at 10:28 UT to at 10:33 UT, as 
seen in continuum intensity. Not only does the diameter vary, 
but also the velocity. This particular UD is interesting because 
it harbors both strong upflows and downflows during part of 
its existence. It starts showing an upflow of -950 m s~^ , but no 
downflow. At 10:31 UT the upflows have increased to -1530 
m s~^ , and a downflow patch of around 700 m s~^ has appeared 
next to the upflow patch. At 10:33 UT the upflows are the 
same while the downflows have increased to 820 m s~^ . At 
10:40 and 10:41 UT, the motions seem to slow down, reaching 
velocities of -950 m s~^ and 180 m s~^ , respectively. The sizes 
of the upflow and downflow patches do not change much in 
this particular example. 

At 10:49 UT, a fast downflow of 950 ms~^ becomes very 
prominent in the lower left part of the image (yellow circles). 



This patch had existed at least during three snapshots, but it 
is in the last one when it shows the highest velocities. The 
downflows have the shape of a collar surrounding a large UD 
with upflows of -780 m s~^ . 

The red circles pinpoint a group of three UDs in the upper 
part of the FOV. They are wefl visible at 10:31 and 10:33 UT 
but have already disappeared at 10:40 UT. At 10:31 UT, the 
rightmost UD shows an upflow of -1280 ms~^ and an adja- 
cent downflow of 470 m s~^ . The other two UDs harbor up- 
flows of -950 ms~^ on average. At 10:33 UT, the downflow 
patch has vanished and the upflows have decreased to -1150 
ms~^ 

Finally, the white circles highlight one of the UDs that con- 
tain a dark lane. The dark lane appears in intensity at 10:31 
UT and disappears at 10:40 UT as the whole structure be- 
comes blurry. In the velocity maps, we can observe narrow, 
elongated upflows that coincide with the position of the dark 
lane. Their velocities are -670 m s~^ at 10:28 UT, -730 m s~^ 
at 10:31 UT, and -860 m s'^ at 10:33 UT. At 10:31 UT, down- 
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Fig. 5. — Stokes /, Q, U, and V profiles emerging from a UD (solid line), its dark lane (thick solid line), and the adjacent umbra (dashed line). The profiles are 
normalized to the quiet Sun continuum intensity. 



flow patches of 560 m appear at the two endpoints of the 
dark lane. One of the patches is as smaU as Of' 14 in diameter. 
Their velocities quickly decrease with time; 2 minutes later, 
at 10:33 UT, the downflows are no longer visible. 

All of these cases corroborate: (1) the ephemeral life of 
the UDs themselves (around 10 minutes), (2) the ephemeral 
life of their substructures (the dark lanes are well defined for 
only around 4 minutes), and (3) the quick change in their ve- 
locities, with up flow/downflow motions lasting only for a few 
minutes. In fact. lSobotka & PuschmannI 12009) report that the 
substructures observed within UDs (i.e., the dark lanes) vary 
with a typical time scale of a bout 3 minutes. According to 
[Sobotka & Puschmann ( 2009), dark lanes disappear and reap- 
pear during the evolution of central UDs. Our time sequence 
is relatively short, and we have been unable to observe the 
reappearance of any dark lane. 

3.4. Polarization signals and magnetic properties 

In Figure m we present the four Stokes profiles of an UD, 
its dark lane, and the surrounding umbra. The selected UD 
is the one enclosed by the white circle in Figure IH With 
a continuum approaching that of the quiet Sun, this UD is 
brighter than the adjacent umbra by a factor of about 2.6. 
[Sobotka et al. (1993) deduced that the brightness of UDs is at 
most 3 times that of the background umbra. Our observations 
appear to confirm this result, derived from a careful analy- 
sis of high-resolution photometric measurements. A factor of 
2.6 is also in agreement with the average continuum ratio of 
2. 58 obtain ed from magnetoconvection simulations at 630 nm 
dBharti et al.ll2010l) . The dark lane, in turn, is observed to be 
about 25% darker than the brightest part of the UD. 

The UD and the dark lane show signifi cantly larger Stokes 
V signals than the umbral background. Socas-Navarro et all 
([2004) pointed out that this does not imply stronger magnetic 
fields, but is just a consequence of a steeper source function 
gradient in hot structures. The UD and the dark lane also show 
enhanced Stokes Q and U signals, suggesting more inclined 



fields. Overall, the spectra dis played in Figu re [5] are very sim- 
ilar to those presented by Soc as-Navarro et al. ( 2004), except 
that these authors did not resolve the contribution of the dark 
lane. 

Figure [6] shows maps of circular and linear polarization 
signals, NCPs, field strengths, and field inclinations in the 
pore. The circles indicate selected UDs, and the dark lanes 
are named DLl and DL2. 

In linear polarization, the UDs stand out clearly above the 
umbral background. The umbra has weak signals due to its 
nearly vertical fields and the small heliocentric angle of the 
observations. UDs 5 and 7 get the higher LP values among 
the selected UDs. Interestingly, the dark lanes DLl and DL2 
exhibit smaller degrees of linear polarization than their cor- 
responding UDs. In circular polarization, the UDs are also 
the most prominent features, despite the fact that the umbral 
background is now very intense. The dark lanes again show 
smaller signals than their UDs. 

It is remarkable that many UDs have strong negative NCPs. 
The temporal evolution of those signals is displayed in the 
lower panels of Figure IH The NCP patches are shifted with 
respect to the photometric center of the UDs, and tend to co- 
incide with downflow areas (e.g., UD4 and UD5) or with the 
border of an upflow patch (e.g., UD6). The dark lanes ex- 
hibit nearly zero NCPs, but the downflows at their endpoints 
sometimes show large negative values (e.g., DL2). The asym- 
metries of Stok es V provide us with valuable information. 
As discussed by Solanki & Montavon ( Il993|) . the sign of the 
NCP is determined by the following rules: 



sign (NCP) = -sign ( Vuuc ^ 



sign (NCP) = -sign V^u^ 



dvLos d|cos7| 
dr dr 



(9) 
(10) 



where Vbiue represents the sign of the blue lobe of Stokes V 
(negative in our case), B is the magnetic field strength (pos- 
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Fig. 6. — Magnetic properties of UDs. In all panels we have marked the edges of the pore with contour lines and a few UDs with circles. Tickmarks 
are separated by I" . Top: continuum image of the pore and its magnetoconvective structures (left); mean linear polarization degree (middle); mean circular 
polarization degree (right). Bottom: net circular polarization (left); magnetic field strength (middle); magnetic field inclination (right). 



itive by definition), and vlos is positive for downflows and 
negative for upflows. We expect larger changes of B with 
height than of 7, hence the first relation above is likely to set 
the NCP sign. The NCP patches coincide with downflows for 
the most part. Since the bisector velocities increase toward the 
continuum, dvLos/dr > in those areas. Thus, to produce a 
negative NCP it is necessary that dB/dr be negative, i.e., the 
magnetic field strength must decrease as the optical depth in- 
creases. This result lends support to the idea that UDs are 
structures with significantly reduced fields near or just below 
the continuum forming layers. 

The interpretation of the observed polarization signals is 
difficult because of the interplay of the various atmospheric 
parameters and the existence of thermal effects. The field 



strength and inclination maps of Figure[6l however, show that 
many UDs are associated with weaker and slightly more in- 
clined fields than their surroundings. The dark lanes have 
even weaker fields. The areas of weaker fields seem to co- 
incide in general with the locations of downflows but not with 
upflows (the dark lanes are an exception, showing both up- 
flows and weaker fields). In the simulations, the downflows of 
the largest UDs are sometimes associ ated with reversed mag- 
netic polarities, but only near r = 1 (iBharti et aLllMToh . We 
do not observe opposite-polarity fields in our observations. 

Table [2I quantifies the properties of the UDs marked in Fig- 
ure [6l To allow comparisons, we also give the properties of 
nearby points in the umbral background (UB). The magnetic 
field strength of the UDs is considerably weaker than that 
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TABLE 2 
Magnetic Properties of Selected UDs 



Structure 


B 


7 


LP 


CP 


NCP 




(G) 


(°) 


(%) 


(%) 


(mA) 


UDl 


1610 


139 


4.0 


9.8 


-1.5 


DLl 


1220 


141 


3.2 


8.1 


-0.5 


UBl 


1900 


150 


2.1 


8.2 


-0.4 


UD2 


1500 


138 


3.8 


9.3 


-3 


DL2 


1020 


141 


2.8 


7.4 


0.8 


UB2 


1800 


142 


2.1 


7.0 


-0.8 


UD4 


800 


124 


4.1 


5.0 


-4 


UB4 


1200 


134 


2.5 


5.0 


-1.5 


UD5 


1190 


129 


4.4 


7.4 


-9 


UB5 


1140 


140 


2.2 


5.0 


-2 


UD6 


1280 


138 


3.7 


6.9 


-4 


UB6 


1450 


141 


2.2 


5.8 


-2.7 


UD7 


1070 


121 


4.9 


6.0 


-3 


UB7 


1480 


140 


2.3 


5.6 


-1.2 



Note. — Column 1: number of UD and adjacent UB in Figure |6] Col- 
umn 2: magnetic field strength. Column 3: inclination of magnetic field 
vector to the LOS. Column 4: mean linear polarization degree. Column 5: 
mean circular polarization degree. Column 6: net circular polarization. 



of the umbra, which peaks at about 2400 G. Differences are 
smaller when the comparison is made with the immediate UB. 
In some cases, like UD4, the magnetic field is as weak as 800 
G. The dark lanes are the regions with weaker fields within the 
UDs: DLl and DL2 show strengths of 1000-1200 G, while 
the bright parts of the corresponding UDs have fields of some 
1500-1600 G. In the umbra, the field is inclined to the LOS 
by up to 160°, as expected for vertical fields at an heliocen- 
tric angle of ~30°. By contrast, the UDs exhibit inclinations 
closer to 90° and thus possess more horizontal fields. UDs 
4 and 7 in particular have the largest inclinations (124° and 
121°, respectively), perhaps because of their location near the 
edge of the pore. Compared with the adjacent background, 
UDs show inclination differences of 3°-17°. The dark lanes 
have similar inclinations as their UDs, but higher than other 
unresolved UDs. 

In summary, UDs exhibit weaker magnetic fields than the 
surrounding umbra as well as more incline d field lines. These 
results confirm pr evious works (e.g. Socas-Navarro et all 
120041 : iRiethmuller et al. 2008; Sobotka & Jurcak 2009). The 
new information here is that of the dark lanes, whose mag- 
netic properties have not been determined earlier. We report 
that the dark lanes exhibit an even weaker magnetic field than 
their associated UDs, up to 500 G less, as well as smaller LP, 
CP, and NCP signals. The field inclination, however, does not 
differ from that of the UDs. Another new result is the exis- 
tence of enhanced NCP in the downflow patches observed at 
the periphery of UDs. The sign of the NCP indicates that the 
field strength decreases with depth in the photosphere. 

4. DISCUSSION AND CONCLUSIONS 

We have presented the first spectropolarimetric measure- 
ments of UDs at a resolution of 0.''14. Our observations re- 
veal the existence of substruct ures within UDs in the form of 
dark lanes. On ly Bharti et al] (l2007bh . iRimmeld (^008), and 
ISobotka & Pus chmann ( 2009|) have detected these structures 
before. Rimmelci (i2008i) estimated their size to be Of' 12 in 
G-band images, righ t at the diffraction limit of the Dunn So- 
lar Telescope, while ISobotka & Puschmann' ("2009) reported 
widths of less than 0."14 from broadband 602 nm filtergrams 
taken at the SST. Our observations also indicate sizes of 0.''14. 



At the resolution of CRISP, however, not all UDs possess a 
dark lane. 

In this paper, we have focused on the velocity field of UDs. 
One can find in the literature a variety of works describing 
the morphological properties of UDs, but very few attempts 
have been made to determine their velocities. However, very 
few attempts have been made to determine their velocities. A 
good knowledge of the flow field of UDs is important both to 
understand the energy transport in strongly magnetized media 
and to discern between models of sunspot structure. 

Measuring UD velocities is a di fficult task, as recognized 
by several authors (e.g.,[R immele 20041: ISocas-Navarro et aP 
120041: ISchUssler & V5gleri>200 6). Our results, and those from 
previous works, point to the existence of large LOS velocities 
and reduced magnetic field strengths, but only in deep lay- 
ers of the photosphere -near the continuum forming region- 
definitely deeper than the layers t raced by the cor e of many 
photospheric lines. For example, IRimmeld ("2004") finds up- 
flows of more than 1 km s~^ in the C I 538.0 nm line (whose 
core forms at approximately 40 km), while he observes veloc- 
ities of less than 300 ms~^ in Fe I 557.6 n m (for mation height 
of around 320 km). ISchussler & Voglerl (l2006h also mention 
the difficulty of observing flows and reduced field strengths 
in UDs because the surfaces of equal optical depth are locally 
elevated. Other factors such as the small sizes and the rapid 
evolution of these structures complicate the determination of 
their velocities and magnetic fields. 

Obse rvational studies ha ve found upfl ows of 100 ms~ ^ 



( Socas- Navarro et al 
i2009» 
ms 



2004. 600 m s:;j_ (|Sobotka & JurcakI 



800 ms-i ('Riethmuller et aL] |2008D, and even 1000 
^ (Rimmele 2004) in peripheral UDs. Our data con- 
firm the existence of relatively strong upflows of around 1000 
ms~\ with peaks up to 1500 ms~^ The upward velocities 
observed in MHD si mulations do not exceed 2000 m s~^ at 
photospheric levels (Bhar ti et aD l2010V 

The most relevant and novel result of this paper is the find- 
ing of downflow s at the edges of some UDs. If the UD has a 
dark lane, then the downflows are observed at its endpoints. 
The downflow velocities range from 400 m s~^ to almost 1000 
m s~^ The existence of downflows seems to fit the scenario 
proposed by Schiissler & Vogler ( 2006), in which UDs are 
the result of magnetoconvection in regions of reduced field 
strengths with both upflows and return downflows. We be- 
lieve this is the first time that the predicted downflows are 
reliably measured within UDs, mostly as a consequence of 
the high resolution provided by CRISP and the SST. We have 
also examined the LOS velocities at different intensity levels. 
Since the 80% intensity level shows the strongest velocities 
and motions at other intensity levels are gradually of lesser 
magnitude, a gradient of velocity with height appears to exist. 
This suggests that the flows occur preferentially in the deep 
layers of the photosphere. 

A comparison between our Figure Of) and Figure 1(b) of 
ISchussler & Vogled (l2006l) reveals some similarities, like the 
fact that the upflows are co-located with the centers of the 
UDs and their dark lanes. Also in good agreement is the 
fact that the downflows are observed around the UD struc- 
ture, but most prominently at the endpoints of the dark lanes. 
However, our measured velocities appear to be slightly larger 
than those predicted by the simulations (see, e.g.. Figure 16 
of Bharti et al. 2010). In the same manner, the temporal evo- 
lution shown in our Figured] may be compared with Figure 5 
of Schiissler & Vogler ( 2006), where the highly dynamic and 
transient nature of the UDs can be seen. 
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We have also determined the magnetic properties of UDs. 
New in this work is the description of the properties of their 
substructure, the dark lanes. UDs exhibit weaker magnetic 
fields and more inclined field lines than the surrounding um- 
bra, confirming previous analyses (e.g. Socas-Navarro et aP 
120041: .Riethmuller et a l. 2008; Sobotka & Jurcak 2009). We 
report a magnetic field weakening of up to 500 G between 
UDs and their adjacent umbral backgro und, similar to the 
values given by iRiethmtiller et aT] (l2008h . According to our 
data, the magnetic field in the UDs is more horizontal than 
in the immediate surroundings by as much as ~20°. Be- 
fore this work, the magnetic properties of the dark lanes were 
largely unknown. We report that they exhibit an even weaker 
field than their associated UDs, up to 500 G less. The dark 
lanes also show smaller linear and circular polarization sig- 
nals. Their inclinations, however, do not differ from those of 
the associated UDs. Another new result is the existence of 
enhanced NCP in UDs. We observe tiny patches of NCP that 
tend to coincide with the downflowing areas. In those regions. 



the NCP sign indicates a reduction of the field strength to- 
wards deeper photospheric layers. This a gain appears to con- 
firm t he results of numerical simulations (ISchtissler & Voglerl 
l2QQ6h . 
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